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Abstract 

Indirect precision data are used to constrain the masses of possible extra Z' bosons 
and their mixings with the ordinary Z. We study a variety of Z' bosons as they appear 
in Eq and left-right unification models, the sequential Z boson, and the example of an 
additional U(l) in a concrete model from heterotic string theory. In all cases the mixings 
are severely constrained (sin# < 0.01). The lower mass limits are generally of the order 
of several hundred GeV and competitive with collider bounds. The exception is the Z^ 
boson, whose vector couplings vanish and whose limits are weaker. The results change 
little when the p parameter is allowed, which corresponds to a completely arbitrary Higgs 
sector. On the other hand, in specific models with minimal Higgs structures the limits 
are generally pushed into the TeV region. 

The possibility of additional neutral gauge bosons, Z's, is among the best motivated types 
of physics beyond the Standard Model (SM). They are predicted by most unifying theories, such 
as Grand Unified Theories (GUTs), left-right unification, superstring theories and their strong 
coupling generalizations. In many cases their masses remain unpredicted and may or may not 
be of the electroweak scale. In the context of superstring models, however, which are much 
more constrained than purely field theoretical models, they are often predicted to arise at the 
electroweak scale, as we will discuss below. In this paper we consider six different types of Z' 
bosons: 

1. The Z x boson is denned by 50(10) — > SU(5) x U(l) x . This boson is also the unique 
solution to the conditions of (i) family universality, (ii) no extra matter other than the 
right-handed neutrino, (iii) absence of gauge and mixed gauge/gravitational anomalies, 
and (iv) orthogonality to the hypercharge generator. In the context of a minimal 5*0(10) 
GUT, conditions (i) and (ii) are satisfied by assumption, while (iii) and (iv) are automatic. 
Relaxing condition (iv) allows other solutions (including the Z^ below) which differ 
from the Z x by a shift proportional to the third component of the right-handed isospin 
generator. 

2. The Z^p boson is defined by E 6 — > 5*0(10) x U{1)^. It possesses only axial-vector couplings 
to the ordinary fermions. As a consequence it is the least constrained of our examples. 
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3. The Zjj boson is the linear combination y3/8Z x — y5/SZ^. It occurs in Calabi-Yau 
compactifications [jl] of the heterotic string if E 6 breaks directly to a rank 5 subgroup [[| 
via the Hosotani mechanism 0]. 

4. The Z LR boson occurs in left-right models with gauge group SU (3)c x SU{2) l X S*?/ (2) fi x 
^(1)b-l C 50(10) and is defined through the current Jlr = ^3/5 [a J 3 r — l/(2a) Jb-l]- 
J 3 r couples to the third component of SU(2) R , B and L coincide with baryon and lepton 

number for the ordinary fermions, a = \jg\jg\ cot 2 0~w ~ 1? where gL,R are the SU(2)l,r 
gauge couplings, and 9w is the weak mixing angle. 

5. The sequential Z$m boson is defined to have the same couplings to fermions as the SM 
Z boson. Such a boson is not expected in the context of gauge theories unless it has 
different couplings to exotic fermions than the ordinary Z. However, it serves as a useful 
reference case when comparing constraints from various sources. It could also play the 
role of an excited state of the ordinary Z in models with extra dimensions at the weak 
scale. 

6. Finally we consider a superstring motivated Z string boson appearing in a specific model || 

based on the free fermionic string construction with real 
fermions. This model has been investigated in consider- 
able detail [0] with the goal of understanding some of the 
characteristics of (weakly coupled) string theories, and of 
contrasting them with the more conventional ideas such 
as GUTs. While this specific model itself is not realis- 
tic (for example it fails to produce an acceptable fermion 
mass spectrum) the predicted Z string it contains is not 
ruled out. Its coupling strength is predicted and so are 
its fermion couplings. It is particularly interesting in that 
its couplings are family non-universal. While this may 
induce problems with too large flavor changing neutral 
currents through a violation of the GIM mechanism 0, 
we will not address this issue here. Another important 
observation is that such a Z string can be naturally at the 
electroweak scale ||. The basic reasons are strong restric- 
tions on the superpotential, and that in a given model 
the sectors of supersymmetry (SUSY) breaking, and its 
mediation are usually not arbitrary. In string models one 

also expects bilinear terms in the Higgs superpotential to vanish at tree level (otherwise 
they should be of the order of the Planck scale) and to be generated in the process of 
radiative symmetry breaking. This is linked to the top quark Yukawa coupling driven 
symmetry breaking and typically involves extra Higgs singlets which are predicted in 
many string models. 

In all cases, there is a relation between the mixing angle 9 between the ordinary Z and the 
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extra Z' from the diagonalization of the neutral vector boson mass matrix 

Ml - Mf = 1 - pp/pi 
M z , 2 - Ml Po /p 2 - 1 



tan 2 ^ = ^S = ^^, (1) 



where Mz and Mz> are the physical boson masses, and Mq is the mass of the ordinary Z in 
the absence of mixing. The second equality in relation ([[]) uses the neutral and charged boson 
mass interdependence which reads at tree level, 

M. = -J^. (2) 



Fp~Z, cos 9\y 

As in the case of the SU(3) C x SU(2) L x U(l) Y model, 

_ -&+**) i<&)i 2 

where ^ (i3j) is (the third component of) the weak isospin of the Higgs field fa. po = 1 if only 
SU(2) Higgs doublets and singlets are present, in which case M would be known independently. 
Nondegenerate SU{2) multiplets of extra fermions and scalars affect the W and Z self-energies 
at the loop level, and therefore contribute to the T parameter |J. They can arise, for example, 
in E 6 models, and in their presence p should be replaced by po/(l — a(M z )T). 

If the Higgs U(l)' quantum numbers are known, as well, there will be an extra constraint, 

where <?i )2 are the U(l) and U(l)' gauge couplings with 

92 = y| sin^VX^i. (5) 

A = 1 if the GUT group breaks directly to SU(3) x SU{2) x £7(1) x £7(1)', while in general A is 
still of 0(1). We will quote our results assuming A = 1, but our limits also apply to y/\ sin 8 
and for other values of A (we always assume M z > 3> Mz). Similar to po, 



' ^itUifaW 



is another function of vacuum expectation values (VEVs), where Q\ are the U(l)' charges. For 
minimal cases, the functions C are given explicitly in Table III of reference fTOfl . Similarly, for 
the fermion couplings of the various Z's we refer to Table II of the same work. 

There is the possibility of an extra gauge invariant term, mixing the field strength tensors 
of the hypercharge and the new gauge bosons. We do not consider such a term here, since it is 
expected to be small in typical modelsfj. The phenomenology of gauge kinetic mixing has been 



x It was shown in Ref. |T^] that a relatively large kinetic mixing term can be generated at the loop level when 
Higgs doublets from a 78 representation of Eq are employed. However, restriction to Higgs doublets from 27 
and 27 representations yields much smaller effects fl3|j. 
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Quantity 


Group (s) 


Value 


Standard Model 


pull 


M z [GeV] 


LEP 


91.1867 ±0.0021 


91.1865 ±0.0021 


0.1 


T z [GeV] 


LEP 


2.4939 ±0.0024 


2.4957 ±0.0017 


-0.8 


r(had) [GeV] 


LEP 


1.7423 ±0.0023 


1.7424 ±0.0016 


— 


r(inv) [MeV] 


LEP 


500.1 ±1.9 


501.6 ±0.2 


— 


r (£+£-) [MeV] 


LEP 


83.90 ±0.10 


83.98 ±0.03 


— 


ow [nb] 


LEP 


41.491 ±0.058 


41.473 ±0.015 


0.3 


Re 


LEP 


20.783 ± 0.052 


20.748 ±0.019 


0.7 


R^ 


LEP 


20.789 ± 0.034 


20.749 ±0.019 


1.2 


R T 


LEP 


20.764 ± 0.045 


20.794 ±0.019 


-0.7 


A FB {e) 


LEP 


0.0153 ±0.0025 


0.0161 ±0.0003 


-0.3 




LEP 


0.0164 ±0.0013 




0.2 


A fb (t) 


LEP 


0.0183 ±0.0017 




1.3 


Rb 


LEP + SLD 


0.21656 ±0.00074 


0.2158 ±0.0002 


1.0 


R c 


LEP + SLD 


0.1735 ±0.0044 


0.1723 ±0.0001 


0.3 


Rs,dJ R{d+u+s) 


OPAL 


0.371 ±0.023 


0.3592 ±0.0001 


0.5 


A FB (b) 


LEP 


0.0990 ±0.0021 


0.1028 ±0.0010 


-1.8 


A FB (c) 


LEP 


0.0709 ± 0.0044 


0.0734 ±0.0008 


-0.6 


A FB {s) 


DELPHI + OPAL 


0.101 ±0.015 


0.1029 ±0.0010 


-0.1 


A b 


SLD 


0.867 ±0.035 


0.9347 ±0.0001 


-1.9 


A c 


SLD 


0.647 ±0.040 


0.6676 ±0.0006 


-0.5 


A s 


SLD 


0.82 ±0.12 


0.9356 ±0.0001 


-1.0 


A LR (hadrons) 


SLD 


0.1510 ±0.0025 


0.1466 ±0.0015 


1.8 


Alr (leptons) 


SLD 


0.1504 ±0.0072 




0.5 


A, 


SLD 


0.120 ±0.019 




-1.4 


A T 


SLD 


0.142 ± 0.019 




-0.2 


A e (Q LR ) 


SLD 


0.162 ±0.043 




0.4 


A T (V T ) 


LEP 


0.1431 ±0.0045 




-0.8 


A e (V T ) 

sXQfb) 


LEP 


0.1479 ±0.0051 




0.3 


LEP 


0.2321 ±0.0010 


0.2316 ±0.0002 


0.5 



Table 1: Z pole precision observables from LEP ||1^|T5| and the SLC fi6|-|I8|] . Shown are 
the experimental results, the SM predictions, and the pulls. The SM errors are from the 
uncertainties in Mz, \nM B , m t , a(Mz), and a s . They have been treated as Gaussian and 
their correlations have been taken into account. The first set of measurements is from the Z 
line shape and leptonic forward-backward asymmetries, A FB (£) = 3/AA e A£. The hadronic, 
invisible, and leptonic decay widths are not independent of the total width, the hadronic peak 
cross section, and the Ri = r(had)/T(£ + £~), and are shown for illustration only. The second 
set represents the quark sector, where R q = r(gg)/r(had), and A q = 4/3A£f(g) is a function 
of the effective weak mixing angle of quark q. The third set is a variety of polarization and 
forward-backward asymmetries sensitive to the leptonic weak mixing angle, where analogous 
definitions apply. For details see reference [|iTfl . 
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reviewed in reference [21]; constraints on leptophobic Z' bosons can be found in reference [22]; 
reference |23] are detailed reviews on extra neutral gauge bosons at colliders; and for related 
topics in neutral current physics see reference |24] . 

After the Z' properties have been specified, the new contributions to the precision ob- 
servables can be computed using the formalism presented in references [|IO| ) p5|] . The new Z' 
boson is treated as a small perturbation to the SM relations. The most important effect is the 
modification of the M^-M^-sin 2 6-w interdependence via Z-Z' mixing. Z pole observables are 
affected by the modification of the Z couplings to fermions, which is another manifestation of 
the Z' admixture. This change in couplings is also relevant for the low energy observables from 
neutrino scattering and atomic parity violation (APV). For these there will be additional effects 
from Z' exchange, and Z-Z' ( r y—Z l ) interference. Such effects can be neglected at the Z pole, 
but the interference terms are relevant for the cross section measurements at LEP 2. They give 
interesting constraints on Mz* [^7| practically independent of the mixing parameter 9. These 
constraints are complementary to the direct search limits at Fermilab p8| , where additional 
assumptions about possible exotic decay channels have to be specified. 



Quantity 


Group (s) 


Value 


Standard Model 


pull 


m t [GeV] 


Tevatron 


173.8 ±5.0 


171.4 ±4.8 


0.5 


M w [GeV] 


Tevatron + UA2 


80.404 ±0.087 


80.362 ±0.023 


0.5 


M w [GeV] 


LEP 


80.37 ±0.09 




0.1 


R~ 


NuTeV 


0.2277 ±0.0021 ±0.0007 


0.2297 ±0.0003 


-0.9 


R v 


CCFR 


0.5820 ±0.0027 ±0.0031 


0.5827 ±0.0005 


-0.2 


R v 


CDHS 


0.3096 ±0.0033 ±0.0028 


0.3089 ±0.0003 


0.2 


R v 


CHARM 


0.3021 ±0.0031 ±0.0026 




-1.7 


RF 


CDHS 


0.384 ±0.016 ±0.007 


0.3859 ±0.0003 


-0.1 


RF 


CHARM 


0.403 ±0.014 ±0.007 




1.1 


R D 


CDHS 1979 


0.365 ±0.015 ±0.007 


0.3813 ±0.0003 


-1.0 


9 v v e 


CHARM II 


-0.035 ±0.017 


-0.0395 ± 0.0004 




9 u v e 


all 


-0.041 ±0.015 




-0.1 


97 


CHARM II 


-0.503 ±0.017 


-0.5063 ±0.0002 




9a 


all 


-0.507 ±0.014 




-0.1 


Qw{Cs) 


Boulder 


-72.41 ±0.25 ±0.80 


-73.10 ±0.04 


0.8 


Qw(T\) 


Oxford + Seattle 


-114.8 ± 1.2 ±3.4 


-116.7 ± 0.1 


0.5 



Table 2: Non-Z pole precision observables from Fermilab [F29|-p3]|, CERN ||14| , p4| , and elsewhere. 
The second error after the experimental value, where given, is theoretical. The SM errors are 
from the inputs as in Table [I]. The various quantities R are cross section ratios from z/-hadron 



scattering, where the CHARM [Q results have been adjusted to CDHS |36| conditions, and 
can be directly compared. The gv,A are effective four- Fermi couplings from v-e scattering ] |37]j . 
and Qw denotes the weak charge appearing in APV [|38]-f4T|]. 



In our analysis we use the data as of ICHEP 98 at Vancouver. It includes the very 
precise Z pole measurements from LEP and the SLC, which are close to being finalized; the W 
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boson and top quark mass measurements, My? and m t , from the Tevatron run I, and further 
My/ determinations from LEP 2; results from deep inelastic i/-hadron scattering at CERN and 
Fermilab; //-electron scattering; and atomic parity violation. The low energy measurements 
in neutrino scattering and APV are very important in the presence of new physics, and in 
particular, for the Z' bosons discussed here. They offer complementary information about Z' 
exchange and interference effects, which are suppressed at the Z pole. The Z pole observables 
are summarized in Table p] and the non-Z pole observables in Table 0. For more details and 
further references we refer to our recent reviews 
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The theoretical evaluation uses the FORTRAN package GAPP dedicated to the 
Global Analysis of Particle Properties. GAPP attempts to gather all available theoretical 
and experimental information from precision measurements in particle physics. It treats all 
relevant SM inputs and new physics parameters as global fit parameters. For clarity and to 
minimize CPU costs it avoids numerical integrations throughout. GAPP is based on the MS 
renormalization scheme which demonstrably avoids large expansion coefficients. 





z x 




Zf\ 


Zlr 


ZsM 


7 

•^string 


SM 






551(591) 


151(162) 


379(433) 


570(609) 


822(924) 


582(618) 






sin 9 


-0.0006 


+0.0004 


-0.0010 


+0.0002 


-0.0015 


-0.0002 






sin 9 min 


-0.0022 


-0.0015 


-0.0058 


-0.0010 


-0.0040 


-0.0011 




Po free 


sin 9 max 


+0.0020 


+0.0021 


+0.0019 


+0.0022 


+0.0008 


+0.0008 




Po 


0.9993 


0.9974 


0.9979 


0.9995 


0.9982 


0.9996 


0.9996 




Po 


0.9931 


0.9923 


0.9931 


0.9917 


0.9933 


0.9986 


0.9985 




Po 


1.0010 


1.0017 


1.0017 


1.0013 


1.0018 


1.0011 


1.0017 




x 2 - 

A mm 


27.62 


27.52 


27.34 


27.71 


26.83 


27.34 


28.37 






545(582) 


146(155) 


365(408) 


564(602) 


809(894) 


578(612) 






sin 9 


-0.0003 


+0.0005 


-0.0026 


+0.0003 


-0.0019 


-0.0002 




Po = 1 


sin 9 min 


-0.0020 


-0.0013 


-0.0062 


-0.0009 


-0.0041 


-0.0011 






sin 6> ma x 


+0.0015 


+0.0024 


+0.0011 


+0.0017 


+0.0003 


+0.0007 






Amin 


28.43 


28.09 


28.17 


28.22 


27.43 


27.82 


28.79 



Table 3: Mass limits [in GeV] on extra Z' bosons and constraints on Z-Z' mixing for two 
classes of Higgs sectors. The upper part of the Table allows p as a free fit parameter and 
corresponds to a completely arbitrary Higgs sector. The lower part assumes p = 1, but is 
arbitrary otherwise. The first (second) numbers correspond to the 95 (90)% CL lower mass 
limits. Below this we show the central values and the 95% lower and upper limits on sin 9. Also 
shown are the central values and 95% limits for po as a fit parameter. Finally we indicate the 
minimal \ 2 for each model. The last column is included for comparison with the standard case 
of only one Z boson. All results assume Mz < Mjj < 1 TeV. 

In Tables [3] and ^ and Figure [I] we present our main results. We list lower limits on Z' 
boson masses for a variety of cases. Note, that the new physics, i.e., the Z's and the extra 
Higgs bosons, decouple and that the SM [Mz> = oo, 9 = 0) is well within the allowed regions of 
Figure [!]. As a consequence a rigorous Bayesian integration over the Mz> probability distribution 
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diverges^. Therefore, we approximate the 95 (90)% CL limits by requiring A% 2 = 3.84 (2.71), 
which is motivated by a reference univariate normal distribution. Similarly, we define the 90% 
allowed region in the (Mz>, sin0)-plane by A% 2 < 4.61, here referring to a bivariate Gaussian. 

The first set of fits in Table |3| is for the most general case of completely arbitrary Higgs 
sector, while the second set is for po = 1. Below the mass limits we show the best fit values and 
95% lower and upper limits on the mixing parameter sin#. Results on po and the x 2 minimum 
are also shown. For comparison we have included the SU(2) x £7(1) case in the last column. 

We note that the Z LR boson is equivalent to the Z x boson with a non-vanishing kinetic 
mixing term, 



smx 



B^Z' 



It can be absorbed by a shift of the gauge couplings proportional to the third component of the 
right-handed isospin generator and a rescaling of the coupling ratio A. In models with given 
Higgs structure the parameter C — > C — J^- sin x is shifted, as well, while po is unaffected. 
The limits on masses and mixings of the Z x and Zlr bosons, shown in Tables [| and |4] (with 
A = 1), are indeed quite similar. 

Our mass limits on extra Z' bosons are somewhat stronger than those from a recent anal- 
ysis [ 43 1 of Z's in supersymmetric Eq models [44j. The differences are due to a slightly different 
and more recent data set in our analysis, different implementations of radiative corrections, dif- 
ferent statistical methods^, and our alternative evaluation of the photonic vacuum polarization 
effects |45| with advantages for global fits. Moreover, reference |43| assumes the SUSY inspired 



range for the Higgs mass, Mh < 150 GeV 





Z x 


Z^p 


Z>q 


Zlr 


(7 = 


1368(1528) 


1181(1275) 


470( 498) 


1673(1799) 


a = 1 


643( 688) 


146( 156) 


1075(1235) 


925( 987) 


a = 5 


1210(1314) 


1393(1581) 


1701(1948) 


980(1076) 


a = oo 


1464(1601) 


1810(2039) 


1985(2277) 


1537(1711) 



Table 4: 95 (90)% CL lower mass limits on specific Z' bosons as they appear in models of 
unification. Assumed are minimal Higgs structures and pq = 1. Note that a is defined differently 
for the Z^p and Z n models, and the Z x and Zlr models, respectively, as explained in the text. In 
particular, the versions of the Z x and Z LR models most often considered correspond to a = 0. 



Table f| lists results for specific Higgs charge assignments as they occur in "minimal" 
models: For the Z^ and Z v bosons, we assume an SU{2) Higgs singlet with a large VEV s to 
ensure Mz> 3> Mz, and in addition a pair of Higgs doublets with quantum numbers as in the 
5 + 5 of SU(5) appearing in the 27 of Eq. They receive VEVs v and v, where the combination 



131 is not well-defined unless a non-trivial 



2 The Bayesian confidence integral in Eq. (4.13) of reference 
Jacobian is implicitly included. 

3 The difference goes beyond the more common choices of Bayesian versus frequentist kind of approaches. 
The authors of reference choose to allow three fit parameters, of which only two are independent. This is 
a problematic procedure when parameter estimation is desired and renders confidence intervals ambiguous. 
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M z [GeV] 

2500 i — . — r 
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1500 - 
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500 
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M z [GeV] 
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2000 - 



1500 



1000 



500 - 




0.005 0.01 -0.01 -0.005 

sin 

M z [GeV] 
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2000 - 



1500 
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500 - 
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o ' — ■ — 1 — 1 — ■ — 1 — ■ — 1 — 1 — ■ — 1 — ■ — 1 — 1 — ■ — 1 — 1 — 1 — 1 — ■ — 1 o 1 — ■ — 1 — 1 — ■ — 1 — ■ — 1 — 1 — ■ — ' 

~0.01 -0.005 0.005 0.01 -0.01 -0.005 

sin 6 




Figure 1: 90% CL contours for various Z 1 models. The solid contour lines use the constraint 
Po = 1 (the cross denotes the best fit location for the po = 1 case), while the long-dashed lines 
are for arbitrary Higgs sectors. Also shown are the additional constraints in minimal Higgs 
scenarios for several VEV ratios as discussed in the text. The lower direct production limits 
from CDF |2?| are also shown. They assume that all exotic decay channels are closed, and have 
to be relaxed by about 100 to 150 GeV when all exotic decays (including channels involving 



superparticles) are kinematically allowed [28 



v 2 + v 2 is fixed by the measured value of the Fermi constant. Thus, these models are described 
by an extra parameter a = \v/v\ 2 which is analogous to tan 2 /? in supersymmetric extensions 
of the standard model. 

The Z x model does not depend on the ratio \v/v\ so that C = 2/yW is predicted. If we 
add another Higgs doublet with quantum numbers like the SM leptons and VEV x, we have 
the extra parameter a = \x\ 2 /(v 2 + v 2 ), and then C E l/\/l0 [— 3, 2]. However, in those models 
of SUSY in which this Higgs doublet is identified with the superpartner of a SM lepton doublet, 
one has to require x = a = to avoid severe problems with charged-current universality. 

The Higgs content of the Z x model can be lifted to an appropriate Higgs structure for a 
Z LR model (LR 1), transforming under SU(2) L x SU(2) R x U(1) B -l as (2, 2, 0) + (2, 1, 1/2) + 
(1,2,1/2). The same definitions and remarks apply, except that here C 6 y3/5 [— 1/a, a]. 
Another possibility (LR 2) is a Higgs sector transforming as (2, 2, 0) + (3, 1, 1) + (1, 3, 1) which 
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results in C = y3/5a. There are no SU(2) triplets in a 27 of Eq, but one might find them in 
string models without an intermediate GUT group, if they are realized at higher Kac-Moody 
levels (k > 1). In the context of SUSY the Higgs fields carrying B — L charge have to be 
supplemented by extra fields with opposite charge to cancel anomalies and (in the case of 
LR 1) to generate fermion masses. 

There is no analog of a minimal Higgs sector for the sequential Z$m- In the Z stT i ng 
model po = 1 is predicted, and even the Z' mass and the Z-Z' mixing can be calculated 
in principle. Universal high scale boundary conditions yield too large a value for 6 and are 
excluded [O . Non-universal boundary conditions on the soft SUSY breaking terms can yield 



acceptable mixing. In any case, the concrete realization of the soft terms depends strongly on 
the SUSY breaking and mediation mechanisms. We parametrize our lack of understanding by 
allowing an arbitrary Higgs sector (except for po = !)■ 

In conclusion, indirect constraints from high precision observables on and off the Z pole 
continue to play an important role for searches for new physics, and in particular for extra gauge 
bosons. The obtained mass limits are competitive with current direct searches at colliders with 
the highest attainable energies. Moreover, no assumptions about the absence of exotic decay 
channels are necessary^. The indirect constraints are even much stronger in specific models 
with known Higgs structure, where lower limits are typically in the 1 to 2 TeV range. Finally, 
Z-Z' mixing effects are severely constrained to the sub per cent level in all cases. 
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